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Abstract 

We present a complete next-to-leading order (NLO) calculation for the total cross section for 
inclusive Higgs pair production via bottom-quark fusion at the CERN Large Hadron Collider 
(LHC) in the minimal super symmetric standard model (MSSM) and the minimal supergravity 
model (mSUGRA). We emphasize the contributions of squark and gluino loops (SQCD) and the 
decoupling properties of our results for heavy squark and gluino masses. The enhanced couplings 
of the h quark to the Higgs bosons in supersymmetric models with large tan j3 yield large NLO 
SQCD corrections in some regions of parameter space. 



I. INTRODUCTION 



In the standard model (SM), only one Higgs doublet is introduced and one neutral Higgs 
boson remains after electroweak symmetry breaking. In the Minimal Supersymmetric Stan- 
dard Model (MSSM) [ij, two Higgs doublets are required to break the electroweak symmetry. 
The two Higgs doublets, 0i and 02, couple to fermions with weak isospin —1/2 and +1/2 
respectively [2|]. After spontaneous symmetry breaking, there remain five physical Higgs 
bosons: a singly charged Higgs boson if^, two neutral CP-even scalars h and H, and a 
neutral CP-odd pseudoscalar A. The Higgs potential is constrained by supersymmetry such 
that all tree-level Higgs boson masses and couplings are determined by just two indepen- 
dent parameters, commonly chosen to be the mass of the CP-odd pseudoscalar {Ma) and 
the ratio of vacuum expectation values of the neutral Higgs fields (tan/? = V2/V1). 

In the standard model, gluon fusion is the dominant process for producing a pair of HiKgs 
bosons via triangle and box diagrams with internal top quarks and bottom quarks [l, 0, UH, 
0, S, 0]. The Yukawa coupling of the Higgs boson to bottom quarks in the standard model is 
proportional to rrih/vsM, where Vsm is the vacuum-expectation value of the Higgs field, and is 
hence very weak. Thus, at the Fermilab Tevatron and CERN Large Hadron Collider (LHC), 
the rate for inclusive Higgs pair production from b quark fusion is small [1, • However, it 
can become significant in supersymmetric models with large tan (3 since the Yukawa couplings 
of the Higgs bosons to the b quark are enhanced by 1/ cos/3. In the Minimal Supersymmetric 
Model (MSSM), the production rate of pairs of the lightest neutral Higgs boson {hh) from bb 
fusion is larger than the rate from the gluon-gluon initial state for tan < 35 and moderate 
values of [M^ ~ 300 GeV), while the rate for pair production of the heavier neutral 
Higgs bosons {HH and AA) from bb fusion is highly suppressed relative to the gluon-gluon 



production mechanism at small tan/9[8|, [llj. This makes it of interest to evaluate higher 
order corrections to the rates. 

One of the most powerful ways to distinguish between the Higgs boson of the standard 
model and those of the MSSM is to measure the trilinear neutral Higg s boson couplings. This 



measurement is extremely challenging in most scenarios jl2|, Il3|, [IJ, Il5|, [IGj . The enhanced 
Higgs pair production rates from bottom quark fusion when tan f3 is large provide discovery 
potential for determining the trilinear Higgs couplings in the MSSM. The rate for Higgs 



boson production from bottom quark fusion [17l| has been cornputed to NNL0[18| and the 
electroweak and SUSY QCD corrections included to one-loop[l9|. Also, the rate for Higgs 
production in association with a b quark has been computed to NLO[20|, |2l|, |22|, |23|, l2J, |25 
including SUSY QCD corrections |26|]. 

In this paper, we present a complete next-to-leading order (NLO) SUSY QCD calculation 
for Higgs pair production via bottom quark fusion. The leading order (LO) process is 
bb (j)(f) (where = h,H,A) and the NLO cross section includes both 0{as) and C(l/A) 
corrections, where A = ln{M^/mb) [23, 0, 29, 3^. The subprocess bg — *• b(f)(f) contributes 



0{1/A) corrections to the leading order cross section for bb (pcj). Theoretical predictions 
depend on the number of b quarks tagged. Here, we consider only inclusive processes in 
which there are no tagged b quarks. The NLO pure QCD corrections have been computed 
in Ref. [lo|. The focus of this work is the inclusion of the 0(a,) SUSY QCD (SQCD) 
corrections, which consist of squark and gluino loop contributions. We also present a detailed 
study of the effects of the SUSY parameters on the production rate in the MSSM and 
mSUGRA models. 

In section II, we review the leading-order cross section for pp 00 from bb fusion. In 
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section III, we provide the complete next-to-leading order (NLO) SQCD corrections from 
gluino and squark loops for bb cjxp production. Numerical results are given in Section IV 
and conclusions are drawn in Section V. In addition, there are three Appendices. Appendix 
A defines the scalar integrals [3l|, l32|, |34| | used in our computation. Appendix B contains 
the coefficients used in computing the virtual SQCD corrections and Appendix C presents 
the Yukawa couplings, Higgs trilinear couplings and Higgs-squark-squark couplings in the 
MSSM. 



II. LOWEST ORDER PRODUCTION IN hh 

The leading order (LO) inclusive cross section for pp — * 00 via bb 
(Tlo = J dxidx2 b{xi)b{x2) + b{xi)b{x2) (rLois,t,u){bb 



IS 



(1) 



where b{x) and b{x) are the LO parton distribution functions for bottom quarks in the proton, 
o"Lo("S,t,w) is the parton level cross section for bb 00 and s,t,u are the Mandelstam 
variables. Fig. [T] shows the tree level Feynman diagrams for bb — 00 in the MSSM. 
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FIG. 1: The lowest order Feynman diagrams in the MSSM for hh 



with (j) = h,H, A. 



We assign momenta to the initial and the final state partons with 

b{pi)b{p2) 0(P3)0(P4) , 



(2) 



and pi + P2 = P3 + Pa- The 066 Yukawa coupling, (yf^^b, and the /i00 and if00 trilinear 
couplings, gh(j,(i, and gH<f>(f), can be found in Appendix C. The lowest order rate can be expressed 
in terms of spinor structures: 



Til = v{p2)Plu{pi) 
TiR = v{p2)Pru{pi) 



(3) 



with Pr l = (1 ± 75)/2. Following the simplified ACOT prescription [35|, l36|, l37|, we take 
TTLi, = everywhere except in the Yukawa couplings. The tree level amplitudes of the s— , t— 
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and u— channels are: 



_ I 9hbb9hcP(p , 9Hbb9H<t>(i, \ ^rj. , rj. ^ c 



= 9lbb\{'"iP'i) hPLu{pi) + V{p2) i)^PRu{pi))5ap 



XL 



M'^J.^ + M'^j.jS^p, (4) 

where a,/? are color indices. The corresponding spin- and color-averaged matrix elements 
squared, including interferences terms, are 



{Re{M^M'J) 



0|2\ _9%b('^_ \ 
6 [t ) 

6 \u V? ) 

0]r40\\ _ 9^bb 



(5) 



6 \ tu J 

11 2 2 
/|^0|2\ _ 9hbb9h<t>4>^ I 9Hbb9H4>4>^ 

ghbbgHbbghH9H4><i> ( {s- Ml) {s - Ml) + MnMnTh^H \ 

3 [is-Ml + iMhThl^-ls-M^ + iMHTHiy^' 

The LO element is then. 

Mo = + M° + M° 
{\Mo\') = (|M,T) + (|M,T) + (|M„Y) + 2(it:e(M,°M0)) . 
Finally, the parton level cross section for inclusive bb — > (fxf) production is 

^LO = J j-^\{\Mo\')d,PS{bb cpcp) , (6) 

where d2PS{bb 00) denotes the integral over the two-body phase space and the factor of 
1/2 is from the identical particles in the final state. 

III. NEXT-TO-LEADING ORDER CORRECTIONS FOR bb 00 

The parton level NLO cross section is 

^■MLoixi,X2,IJ,) = ^L0ixi,X2,ll)+SaNL0ixi,X2,IJ,) 

= ^Loixi,X2,iJ,) 1 + Sqcb{xi,X2,ij,)+Ssqcb{xi,X2,ij,) 
Sa^Lo{xi,X2,ii) = S&qcd{xi, X2, IJ,) + S&sqcd{xi,X2, IJ,) , (7) 
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where a-Lo{xi,X2, /i) is the leading order (Born) cross section and Sa-^i,o{xi,X2, /i) is the next- 
to-leading order correction to the Born cross section, Xi^2 are the momentum fractions of 
the partons and fj, = fipt is the renormalization scale. The NLO correction 6aNLo{xi,X2, fJ^) 
contains both the gluon QCD correction, Saqcui^i, X2, and the gluino-squark SQCD 
correction, SasQCD{xi,X2, fi). The gluon QCD correction includes the 0{as) corrections, 
which contain virtual and real gluon emission contributions, as well as the 0(1/ A) corrections 
from the bg hipcf) subprocess. The gluino-squark SQCD correction contains 0{as) gluino- 
sbottom loop contributions. 

Unlike down-type quarks, which only couple to the down-type Higgs at tree level, the 
down-type squarks also couple to the up-type Higgs boson. This leads to mixing in the 
sbottom mass matrix in the hi-, basis. 



m 



m 



bi,2 



bL 



mb{Ah — yUtan/3) 



■miy{Af, — yUtan/3) 



m 



(8) 



where 61^2 are the mass eigenstates. The mass eigenstates m^^ and rrii^ (m^;^ < mi^) are 
defined in terms of a mixing angle. 



bi 



cos 61 sin 01 \ I h 
-sin 61 cos Of, [Ir 



(9) 



The one-loop virtual gluino-sbottom Feynman diagrams for bb —>■ (fxp are shown in Fig. [2l 
The amplitudes corresponding to each diagram in Fig. [2] are computed analytically and 
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FIG. 2: One-loop SQCD virtual corrections to bb — > (p(p. 
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all tensor integrals are reduced to linear combinations of one-loop scalar functions. We 
sum all of the virtual amplitudes to obtain the contribution to the matrix element from 
gluino-sbottom loops in terms of the standard matrix elements of Eqs. [3] and HJ 



rO 



\i=l,3,5,7 1=1,3,5,7 i=2,4,6,8 

10 10 \ 

+ J2 ^i2MuR + E WlL + E ^^4^1^ , (10) 



i=2,4,6,i 



i=l 



2 = 1 



where i corresponds to the numbering of the diagrams in Fig. [2l The analytic expressions 
for each coefficient X can be found in appendix B. For simplification, we define: 



XtL 
Xr 
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1=1,3,5,7 
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1=2,4,6,8 
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E^^3 

i=l 



XtR '- 
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=1,3,5,7 
E 



X 



1=2,4,6, 
10 

i=l 



ii ■ 



(11) 



The one-loop SQCD matrix element is then 



[12) 



Only ultraviolet (UV) divergences occur in the SUSY QCD corrections from the massive 
gluino and squark loops. These UV divergences are removed by the renormalization of the 
bottom quark wavefunction and propagator and the renormalization of the bottom quark 
mass in the Yukawa coupling 26, IsH, [s^, 0, 41, 42|. We use the on- shell renormalization 
scheme of Ref. 12] and define the h quark self energy as. 



S^(p) = #(sUp')-S^(p')75)+m,S|(p2) 

which gives the renormalized propagator, 

Ken = {i>-m,) (t" + + m J S 



6mb\ 



rrib 



(13) 



(14) 



The on-shell renormalization condition requires, 

^ = mfe) = 



lim 



0. 



nib 



(15) 



Computing the bottom quark self energy from gluino-squark loops and ignoring contri 
butions suppressed by powers of nib, we find: 

"Stt 

i f ^ 1 sin 



5i(p, m-g, mi J + Bi{p, m-g, m^J 



Stt \ mb 



Bo{p, m-g, m^J - Bo{p, m-g, m^J 



(16) 
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which yields, 



Smi, 
rrib 



V \p'^=m'^ 



(17) 



The couphng of the b squark to the up-type Higgs doublet induces a modification of the 
tree-level relation between the bottom quark mass and its Yukawa coupling. At large tan /3, 
we absorb this modification by redefining the bottom quark mass occuring in the Yukawa 
couphng, [ii, m, S, S, 113 

nib 

nib 



where 



A. 



— nig^itaxi[jl{nii mi m-g) . 



and the auxiliary function is defined as, 

1 



Jfa, 6, c) 



„2 U-2 2 

-,212 1^ " , 1,2^2 i„ " , ^2^2 ^ ^ 



(n^h^ In — + In — + c^o^ In — ) 



(18) 



(19) 



(20) 



The contributions to the bottom Yukawa couplings which are enhanced at large tan /3 can 
be included to all orders by making the following replacements (lol . 43 . 44 . 45 . 



Qhbb 
QAbb 



Qhbb- 



9Hbb 



L — Afc/(tan/3tana) 

l + Ab 
1 + Afotana/ tan/3 



gAbb- 



l + Afe 
Afo/tan^/? 
1 + A. 



We use these effective Yukawa couplings in our NLO computation. To avoid double-counting, 
we add additional mass counterterms. 



6ni^ 
nib 



aJi + 



tan a tan P 



6ni^ 
rrib 



Ab 



tana 
tan l3 



6m^ 
rrib 



Ab 



(21) 

for h, H and A production, respectively. 

The counterterms arise from the b quark wavefunction renormalization (a factor of SZb/2 
for each external b quark), the b mass renormalization in the (pbb couplings (this occurs twice 
in the t— and u— channel diagrams, and once in the s— channel diagram), and a factor of 
—SZv from the mass renormalization on the internal b quark propagators shown in Fig. [S] 



M, 



CT 



^)(M° + M°) + (5Zy + 
nib / V 



Snib 5m 



mb mb / ' " ■ V mb 

The complete one-loop SUSY QCD contribution is, 

^SQCD ^ Mo + M^^^^ + McT , 



+ 



rrib 



(22) 



(23) 
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FIG. 3: t— and u— channel bottom quark propagator counterterm diagrams. 

and is of course finite. Finally, the SQCD contribution to the matrix element squared is, 
|^SQCD|2 ^ I |2 +2i?e(Mf^^^Mo) + 2Re{McTM)) 

iil'b iil'b 

6mb 6mt s 



+2{5Zv + + "-^)\Ml 
nib rrib 



0|2 



(24) 



with the spin and color averaged result 

47ra;, 



2{Re{M^^^^ Mo)) 
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+ {XtR + XuR + XtL + XuL)ReiM?Ml) 
+ 2Re[{XL + XR)X^*]s} . 



(25) 



IV. RESULTS FOR HIGGS PAIR PRODUCTION FROM BOTTOM QUARK FU- 
SION 

In this section, we present the next-to-leading-order inclusive cross sections for the pro- 
duction of a pair of neutral Higgs bosons via bottom quark fusion in the MSSM and 
mSUGRA models at the CERN LHC. As in our previous paper 13], we use the lowest 
order CTEQ6L1 parton distribution functions (PDFs) [47] with a factorization scale fip and 
the leading-order evolution of the strong coupling as{fiR) with a renormalization scale fiR to 
calculate the LO cross sections and use the CTEQ6M PDFs with the next-to-leading-order 
evolution of asifJ^R) to evaluate the NLO inclusive cross sections. For simplification, we use 
the same renormalization and factorization scales fip = fiR. We evaluate the bottom quark 
mass occuring in the (j)bb Yukawa couplings using the MS mass, rribifi), with a two- loop 
heavy quark running mass with mb{pole) =4.7 GeV and the NLO evolution of the strong 



50|. The Higgs 



coupling constant, modified to decouple the effects of the top quark[48|, |4t ^ 

couplings are in Appendix C 0, 51, 52, 53] and we compute the Higgs boson masses to 
one-loop accuracy |54| . 



A. Results in the Minimal Supersymmetric Model (MSSM) 

In Fig. m we show the LO and NLO cross sections versus the pseudoscalar Higgs mass 
Ma- We assume rrig = rrif^ = m^-^ = —A^ = = Msusy and compute the h squark masses 
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FIG. 4: Next-to-leading order cross sections aNLo{pp 4>4> + X) (solid) and leading order cross 
sections (JLoiPP 4>4> + X) (dash) for Higgs pair production from bottom quark fusion versus the 
pseudoscalar Higgs mass {Ma) with \/5 = 14 TeV and IJ-R = IJ'F = The NLO cross sections 

include both the pure QCD and the SQCD corrections. We use tan /? = 50 (red) and tan /? = 10 
(blue) with MgusY = 1000 GeV. Also shown are the cross sections of gg fusion (dot). 

and mixing angles from Eqa. [8]and[9l Our NLO cross sections include the 0{as) corrections 
from the hb initial state and the 0(1/A) corrections from the hg initial state, along with 
the SQCD corrections from gluino-sbottom loops. We show our results with tan/? = 50 
(red) and tan/5 = 10 (blue) at Msusy = 1000 GeV. To compare with gg fusion, we also 
plot the cross section from the gg initial state (dot) js], @]. We note that the cross sections 
for Higgs pair production in the MSSM are significantly larger than in the standard model, 
due to enhancements at large tan/3. The resonant enhancements due to s— channel scalar 
exchange are clearly visible in the hb hh and bb HH curves. At tan/? = 10, the cross 
section for pair of the lightest neutral Higgs boson {hh) from hb fusion is much larger than 
the cross section from the gluon-gluon initial state, while the rate for pair production of the 
heavier neutral Higgs bosons {HH and AA) from hb fusion is highly suppressed relative to 
the gluon-gluon production. But at tan (3 = 50, gluon fusion dominates pair of the lightest 
neutral Higgs boson {hh) production and is comparible with bb production for the heavier 
neutral Higgs bosons {HH and AA) production. 

Figs. [5] and Fig. [6] show the NLO cross sections versus Ma with Msusy = 1000 GeV and 
tan (3 = 50. We present the NLO cross section with only gluon QCD corrections (dash-dot- 
dot, green), NLO cross section with only gluino-sbottom SQCD corrections (dot, blue), and 
the complete NLO cross section with QCD and SQCD corrections together (dash, red). We 
note: 

• The NLO SQCD correction to bb hh is small. The dominant contribution to the 
NLO correction to light Higgs pair production is from the pure QCD contribution. 

• For heavy Higgs pair and pseudoscalar Higgs pair production, the SUSY QCD cor- 
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FIG. 5: Next-to-leading order cross section aNLo{pp hh + X) in bottom quark fusion versus Ma 
with = 14 TeV, tan /? = 50 and MsusY = 1000 GeV. We plot the NLO cross section with only 
gluon QCD correction (dash-dot-dot, green), NLO cross section with only gluino-sbottom SQCD 
correction (dot, blue) and complete NLO cross section with both gluon QCD and gluino-sbottom 
SQCD corrections (dash, red). 



rections become dominant. The pure gluon NLO contribution is much smaller in 
magnitude than the contribution from SQCD. 

In Fig. [TJ we plot the ratio of the NLO SQCD correction normalized to the Born cross 
section, 5sqcd, with Ma = rrig = —A^ = 1000 GeV, and tan/3 = 10 (dash-dot-dot, green), 
tanP = 30 (dash, blue) and tan/3 = 50 (solid, red). In the limit of large squark masses, 
SsQCD approaches a common non-zero constant for HH and AA production. Light Higgs 
pair production, however, decouples for large M4 and large SUSY masses. This decoupling 



behaviour is also seen in the decay h ^ hh [45|, 155[ and the production process hg —>■ 6/i[2 



In Fig. [HI we fix Ma and all squark masses to be 1000 GeV and plot Ssqcd versus the 
gluino mass nig with tanP = 10 (dash-dot-dot, green), tanP = 30 (dash, blue) and tanP = 50 
(solid, red). This figure does not demonstrate a decoupling behaviour. 



B. Results in the Minimal Supergravity Model (mSUGRA) 

In this model, supersymmetry is assumed to be broken in a hidden sector consisting 
of fields that interact with the usual particles and their superpartners only via gravity. 
Supersymmetry breaking is communicated to the visible sector via gravitational interactions. 
Within the mSUGRA framework, it is assumed that at some high scale (frequently taken to 
be ~ Mgut) all scalar fields have a common SUSY breaking mass Mq, all gauginos have a 



10 



tanp = 10, M 



SUSY ■ 



1000 GeV 



0.15 



X 

+ 

-©- 
■o- 



0.1 



0.05 



il I I I I I I I I I I I I I I I I I 

ha A A ■ 

.\\ Complete _ _ . 

_i\ gluino/squark only 

. ■ \ tiluon only 

■ \\\ 

■ \\\ 
. 

\\ 

~ \\ 
\\ 
\\ 



I I I 



I I I I I I I I I I I I I I I I I I 

HH ■ 



Complete 
gluino/squark only 



gluon only 



_/ ^ 
/ \ 

1 ¥ 

i- \^ 
1 v 



I I I I I I I I I I I r 



'T"kH--i 



150 200 250 
M^(GeV) 



300 150 200 250 
M^(GeV) 



300 



FIG. 6: Next-to-leading order cross sections for (a) (Tnlo{pP HH+X) and (b) (7nlo{pP — * AA+ 
X) from bottom quark fusion versus Ma with ^fS = U TeV, tan /? = 50 and MsusY = 1000 GeV. 
We plot the NLO cross section with only gluon QCD correction (dash-dot-dot, green), NLO cross 
section with only gluino-sbottom SQCD correction (dot, blue) and complete NLO cross section 
with both gluon QCD and gluino-sbottom SQCD corrections (dash, red). 
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FIG. 7: SsQCD versus squark mass rriq with Ma = rrig = 1000 GeV and = 14 TeV for (a) 
bb —>■ hh, (b) bb —>■ HH and (c) bb —>■ AA. tan /? = 10 (dash-dot-dot, green), tan /? = 30 (dash, blue) 
and tan/? = 50 (solid, red). 



common mass M1/2, and all soft SUSY breaking scalar trilinear couplings have a common 
value Aq. Electroweak symmetry breaking is assumed to occur radiatively. This fixes the 
magnitude of the superpotential parameter /x. The soft SUSY breaking bilinear Higgs boson 
mass parameter can be eliminated in favour of tan /3, so that the model is completely specified 
by the parameter set: 

Mo, Ml/2, ^0, tan/3, sign{fi). (26) 
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FIG. 8: 5sQCD versus gluino mass, m^, with Ma = f^b£ ^ ~ ^ ~ 1000 GeV and = 14 TeV for 
(a) bh hh, (b) bb — > HH and (c) bb — > AA. tan/3 = 10 (dash-dot-dot, green), tan /? = 30 (dash, 
blue) and tan/3 = 50 (solid, red). 
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FIG. 9: onlo versus Mo with ^/S = IA TeV and tan (5 = 10. Three cases are shown with /x > 
for Ml/2 = 200 (blue)and 400 GeV (red), /i < for M1/2 = 200 GeV(green). Also shown are the 
LO cross sections, cjlo (dot). 



All the sparticle masses and couplings required for phenomenological analysis can be ob- 
tained via renormalization group evolution between the scale of grand unification and the 
weak scale. 

We show the LO (dot) and NLO (solid) cross sections versus Mq with = for tan (3 = 
10 in Fig. [9] and tan /3 = 50 in Fig. [TOl We plot six curves in each frame, three LO cross 
sections (dot) and three NLO cross sections (solid). The NLO curves include both the pure 
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tanp = 50, A. = 
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FIG. 10: (Jnlo versus Mo with \/5 = 14 TeV and tan (3 = 50. Three cases are shown with /i > 
for Ml/2 = 200 (blue) and 400 GeV (red), // < for M1/2 = 200 GeV(green). Also shown are the 
LO cross sections, cjlo (dot). 

QCD and the SQCD contributions. As shown in the graphs, blue lines are cross sections 
with Ml/2 = 200 GeV, ji < green lines have M1/2 = 200 GeV, /z > 0, and red lines are for 
Ml/2 = 400 GeV, fi < 0. Comparing Figs. [9] and [10] we see a strong dependence on tan/3. 

• When tan (3 = 10, flipping the sign of /i has little effect on either the LO or the NLO 
cross sections with M1/2 = 200 GeV. As tan/3 increases to 50, we notice that flipping 
the sign of fj, has a large effect when Mq < 1200 GeV as shown in Fig. [TOl 

• For large Mq, the cross sections approach a common value, independent of M1/2 and 
II. 

V. CONCLUSIONS 

In this paper we presented the complete 0{as) SUSY QCD corrections to neutral Higgs 
pair production from bottom quark fusion in the Minimal Supersymmetric Model and the 
Minimal Supergravity Model at the CERN LHC. The effects of the SUSY QCD corrections 
from sbottom and gluino loops are significant in some regions of parameter space and the 
decoupling of the SQCD effects is only recovered for light Higgs pair production in the large 
Ma limit. 
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APPENDIX A: THE SCALAR INTEGRALS 



IQ^2 uv/^, 1, J ^2n)^ {k^ - ml)[{k + pY - ml] 



'■2J 

^ ^ . . f d-k 1 

J {2n)'' [k^ - mi)[{k + piy - m2\[[k + pi + p2y - mi] 



167r2 



' ACii{pi,P2,mi,m2,m3)+ hCi2{Pi,P2,mi,m2,m3) 



167r2 

d^'k ]k 



J f27r)" 



(27r)" - mf ) [(A; + - mi] [(A; + pi + ^2)^ - mj] 



A-Dll(pi, P2,P3, 7712, ^3, m4)+ i'2-Di2(pi , P2, Ps, "^i, 7712, ^3, 77^4) 



167r2 

+ hDi3{pi,P2,mi,m2,m3,m4) 



(27r)'* (/c2 - mf ) [(/c + pxf - ml] [{k + pi + ^2)^ - mi] [(A; + pi + P2 + Pa)^ - ml] 
-D123 (pi , P2 , P3 , m-i , 7712 , m3 , 7714) 

= Dl2(pi,P2,P3, mi, 77^2, 7713,7714) - -Dl3(pl,P2,P3,mi,77l2,m3,77l4) 



APPENDIX B: X COEFFICIENTS 

1. X coefficients for hh hh and bb HH 



X12 

Xi3 



87r2 
1 



-Bi(-pi +P3,mp,77igjsin2 5)- _^ _^ ^3 cos^ 6lg 

Bi{-pi +P3,mg, mi^) cos^ + Bi{-pi + P3, m^, 777^ J sin^ 9-^ 



1 S'#6 



^ ma 

^13 



-Bo(-Pi +P3,mg,mi^) - Bo{-pi +p3,mg,m 



62 ' 



sin ^^i; cos 9^^ 



X2i = ^ ii,P3 ^ P4) i = 1, 2, 3, 4 
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^31 = ■^^^mgCo{-pi,p3,mg,mi mr)sm9^,cos9^, 
+ ^nigCo {-Pi,P3,mg,mi mi ) cos^ 9^ 

OT^ g<l>bb 

- ^ m-gCoi-pu Ps^m-g, mi mi ) sm Oi 

- ^ ^m-gCo {-pi,p3,mg,mi mi ) sin $1 cos 9i 

ovr g^bb 

^32 = / ^'^^^ mgCo(-pi, p3, mg, mi , mi ) sin 9i cos 

oTT g^bb 

- rrigCoi-pi, P3,mg, mi mi ) sm 9i 

OTT g4>bb 

+ rngCo{-pi,P3,mg,mi^,miJcos 9i 

oT^ gcfebb 

- ■J-^^^^Co{-pi,p3, mg, mi mi ) sin 9i cos 9i 
oTT gt(,bb 

^33 = - -^9^bbg<j>iiCi2{-Pi, P3, rrig, mi^, mij cos^ 9i 
5#65<^i2 Ci2 {-Pi,P3,mg,mi^, mi^ ) sin cos 
^ 9<l>bb9^i2Ci2 i-Pi,P3,mg,mi^, m^^ ) sin cos 

- ^ ^066^<^22 Ci2 ( -pi ,P3,mg, mi^ , mi^ ) sin^ 9i 
^34 = -^c/#6^0iiC'i2(-pi,P3, mg, mi^,mi^) sin^ 

- ^ 9<i>bb9ci>i2Ci2 {-Pi,P3,mg,mi^, mi^ ) sin cos ^5 

- ^ 9^bb9^i2Ci2 hPi,P3,mg,mi^,miJ sin cos 

- ^ 9<i>bb9 ^22 Ci2{-Pi,P3,mg,mi^, m^^ ) cos^ 9i 

= X3j(t -s-^- ii,p3 ^ P4) i = 1, 2, 3, 4 
1 

^51 = -^-^rngCo{p2,-P4,'mg,mi mi)sm9icos9i 
-o(p2, -P4, mg, mi^,mi^) sin^ 
-o(P2, -P4, mi^,mij cos^ 
Co(P2, -P4, mg, mg^, mgj sin 6lg cos 6lg 



+ 
+ 



1 


9<j)ii 


87r2 


9<pbb 


1 


fl'012 


87r2 


5'</)66 


1 


9^12 


87r2 


9cl)bb 


1 


9(f>22 


87r2 gj,bb 
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^52 = -^^mgCo(P2,-P4,mg,mr ,mr )sin6'^cos6'^ 
+ V^^^9^oiP2, -P4, m^, ) cos^ 

1 9<l>12 ^ , \ ■ 2a 

- ^mgCo{p2, -P4, rUg, mi , ) sin cos 9^ 

OTT^ gcl,bb 

^53 = --^g^bbg4>llCl2iP2,-P4,mg, 7711^,7711^8111^61 

- g<t>bbg^i2Ci2 ip2 ,-p4,mg,mi^, ) sin Oi cos ^5 

- ^ g<t>bbg ,pi2Ci2 {P2 ,-PA,mg, 7711^ , 7711^ ) sin COS 

- gci>bbg<t,22 C12 (P2 , -P4 , m~g , mi^ , mi^ ) cos^ 9^ 
^54 = -^5'#65'</.iiCi2(p2, -P4, m-g, mi^,77iij cos^ 6*^ 

+ ^ 50665<^i2 C12 (p2 , -P4 , , mg^ , mg^ ) sin cos 9i 
+ ^ g<t>bbg^i2Ci2 {P2 , -Pa , m-g , 7711^ , 7711^ ) sin 9^ cos 9^ 

- g<t>bbg ^22 C12 (P2 , -P4 , m~g , mi^ , mi^ ) sin^ 9i 

X(ii = Xzi{t <r^u,p3^ P4) i = 1, 2, 3, 4 

1 q"^ 

X71 = — ^tDi23(-pi,P3,P4,mp,mg^,m^^,m^Jcos^e^ 

- ^ ^^i^^i-Di23 {-Pi,P3,P4,mg,mi^, mg^ , m^^ ) sin cos 

9(pbb 

+ ^ ^123 {-Pi,P3,P4,mg,mi^, mi^ , ) cos^ 9i 

_ 1 g^22g^i2 _^jj^^^ ^ ,P3,PA,mg, mi^ , mi^ , m^^ ) sin 9^ cos 9i 
9^bb 

- 2 ^^123 ( -Pi , P3, P4 , , , m^^ , m^^ ) sm cos 9i 

9(t>bb 

+ ^ ^^-Di23 {-Pi,P3,P4,mg,7ni^,mi^,7niJ sin^ 

o'^ 9<pbb 

1 9ij)22g<t>i2,j^ I \ ■ a a 

- ^ 2 ^^123 ( -P\ , P3, P4 , , m^^ , , m^^ ) sm 9^ cos 

°^ 94,bb 

+ ^*-Di23 {-Pi,P3,P4,mg,7ni ) sin^ 
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+ 
+ 
+ 

+ 
+ 

+ 



o"^ 9ct>bb 

+ ^ ^^i^^i-Di23 {-Pi,P3,P4,mg, mi_^ , r/i^^ , m^^ ) sin $1 cos 9^ 
-Y^tDi23 {-pi,P3,P4,mg,mi^,mi^, m^^ ) sin^ 9^ 

2 ^-^123 ( -Pi , P3 , P4 , mg , , mi^ , mg^ ) sm 9i cos 

2 ^-^123 {-Pi:P3:P4:m~g,mi^, mg^ , m^^ ) sm 9i cos 

1 0'^ 

— 4^t-Di23 {-Pi,P3,P4,'mg,mi^,mi^, m^^ ) cos^ 
2 ^^123 ( -Pi , P3, P4 , , , , m^^ ) sm 9^ cos 

— -^tDi23 {-pi,P3,P4,mg,mi^,mi^, m^^ ) cos 

^73 = -^^Jiim^L'o(-pi,P3,P4,mg,m^^,m^^,m^Jsin^^cos^^ 

- ■^9<i>ii94>i2mgDo{-pi,ps, p^, rrig, mg^ , m^^ , mgj cos^ 9^ 

- -^9li2mgDo{-Pi,P3,P4,mg,mi^,mf^^,miJ sin cos 9^ 
■^9<i>229^i2mgDo{-pi,P3,P4, rUg, m^^, m^^, m^J cos^ 9^ 
■^9<t>i^9^i2mgDo{-pi,ps,p4, rrig, mi^,mi^,mi^) sin^ 
■^9li2'rngDo{-pi,p3,P4,mg,mi^,mi^, m^^ ) sin cos 
■^9<f229c^i2mgDQ{-pi,p3,p4, rUg, mi^,mi^,mij sin^ 
^ ^522^9-^0 i-Pi,P3,P4,mg,mi^, , m^^ ) sin cos 

^74 = --^gliimgDo{-pi,p3,p^,mg,mi^,rni^,miJsm9icos9i 
■^90ii9<pi2mgDo{-pi , p3, P4,mg, m^^ , m^^ , m^^ ) sin^ 
^4i2^fl^o(-Pi, P3, P4, m^^ , m^^ , ) sin cos 9^ 
■^9<p2294>i2mgDo{-pi,p3,p4, rUg, m^^, m^^, m^J sin^ 9i 
■^9<l>ii9cl>i2mgDQ{-pi,ps,p4, rrig, mg^, m^^, mgj cos^ 9i 
9li2^gDo {-pi,P3,PA,mg,mi^, m^^ , m^^ ) sin 9i cos 



+ 
+ 
+ 
+ 



+ 



+ 



'b 

'b 
\ 
'b 

h 
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87r2 



9<p2294>i2mgDo{-pi ,P3,P4,mg, m^^ , m^^ , m^^ ) cos^ 6^ 



+ ^ 422^5-^0 {-Pi,P3,P4,mg,mi^, , m^^ ) sin 6*5 cos 6*^ 

= Xri{t ^ u,p3'^ P4) i = 1, 2, 3, 4 
1 rUn 



X, 



93 = - ^QhughH (g _ M^) + iMfoFft ^° ^ ' ^ ' ' ^ ^ 

9hi2gh4>4>T- ,^2\ I ^- 71^ r 



+ ^ghi2gh<pc^-i—jpr^—j^Co{-pi,pi +P2, rrig, m^^, J sin" 



1 



+ ■^gh22ghH (^^ _ j^2^ ^ iMhTh '^°^~^^'^^ + P2:m~g,mi^,m~Jsm9iCOs9i 



1 77T'"' 

+ ^j7M2gfe<A^ ^^ _ ^2^^ iMfeF;, ^"^"^''^' +P2,mg,m^^,m^Jsin2^^ 

X 777-" 

- ^9hi2ghHJ^zrj^Tf^~^^^Co{-Pi,Pi +P2,m~g,mi^,m~Jcos''e~^ 

1 TTT-" 

+ ^9h22gh<i>^j^^^j^^j^^-^^Co{-pi,pi + p2,mg,mi^,m~Jsmeicosei 

^lOi — ^9i{ghjk ^ gHjk, ghcfxp ^ gH<j><j>) ^ = 3, 4; j,k — 1,2 



2. X coefficients for bb — > AA 



X12 
Xiz 



87r2 
1 



Bi{-pi +P3,m~g, mgj sin^ 6'g + Bi{-pi +ps,m~g, mi J cos^ 9^ 
Bi{-pi + P3, "^5, m^J cos^ 6*^ + + p3, mg, m^J sin^ 9^ 



87r2 t 

Xi3 



Bo{-pi +P3:m~g,miJ - Bo{-pi+p3,m~g,mi^) 



sin 9i cos 9i 



X2i = <r^u,ps'^ P4) i = 1, 2, 3, 4 



^ 1 gA12 ^ I \ 2 a 

^31 = ^ mgCo{-pi,P3,mg,mi mijcos 9^ 
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oT^ 9 Abb 

^ 1 9a\2 ^ , \ ■ 2 a 

^32 = ^ mgCo(-pi,P3,mg,m5 Jsm 

, 1 9 AVI ^ , \ 2n 

^33 = -^fi'Afe65'A12Ci2(-pi,]?3,?"9'"^6i'"^62)^i^^fe'^°^^fo 

+ ^5'A665'Ai2Ci2(-pi, P3, "^p, m^^, m^J sin 61^ cos 61^ 

^34 = --^gAbb9Ai2Ci2{-pi,P2,mg,mi^,mi^)iii\ieicosei 

+ ^5'A665'Ai2Ci2(-pi, P3, "^p, m^^, mgj sin 6lg cos 6lg 

= X^i{t <r^u,p3'^ P4) i = 1, 2, 3, 4 



1 


fi'A12 


87r2 




1 


5'A12 


Stt^ 


9Abb 


1 


9A12 


87r2 


9Abb 


1 


9A12 



+ 

X52 = ■^'^^mgCo{p2,-PA,mg,mi^,miJcos^ei 



+ ^ mgCo{p2:-P4:mg,mi^,mi^)sm Oi 

oT^ 9 Abb 

^53 = -^9Ahh9Ai2Ci2{p2-, -p4, m-g, m^^, m^J sin 6*^ cos 6'^ 

- ■^9Abb9Ai2Ci2{p2, -pA,m~g,mi^,mi^) sin cos 
^54 = -^9Abb9Ai2Ci2{p2, -P4, mp, mi^,mij sin 6*^ cos 

- ■^9Abb9Ai2Ci2{p2, -P4, m-g, mi^,mi^) sin 6*^ cos 6*^ 

Xei = X5j(i it,p3 -s-^- P4) i = 1, 2, 3, 4 



1 


fl'ii2 


87r2 


9Abb 


1 


9ai2 


87r2 


9Abb 


1 


9ai2 


87r2 


9Abb 


1 


9ai2 



+ ^ -;^tDi23 ( -Pi , P3 , P4 , , m^^ , m^^ , m^^ ) sin^ 
^72 = ■^^^tDi23{-pi,p3,p4,mg,mi^,mi^,miJsm'^ei 
+ ■^'^tDi23{-pi,p3,p4, rrig, mi^, m^^, m^J cos^ Oi 

°^ 9 Abb 
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^73 = --^9Ai2'n^9Do{-Pi,P3,P4,rng,mi^,mi^,miJsmeicosei 



87r2 

1 .2 



+ ^^Ai2^fl-Oo(-Pi, P3, P4, mp, m^^, mg^ , m^J sin 6*5 cos 6^ 
= -^c/Ai2^9-Co(-pi,P3,P4,mg,m5^,m^^,m5jsine^cose^ 



87r2 

+ Ai2^fl-Do(-pi, P3, P4, mg, mg^, , m^J sin cos 6^ 



87r2 

Xgj = X'ri{t <r^u,ps-^ P4) i = 1, 2, 3, 4 

^93 = -^^fciiffftAA ^^ _ ^ iMhTh ^^^~^^'^^ ^bi^ si^^bcos^b 

- ^gfei2gfe^A ^^ - M^) + iMhVh ^''^~^^'^^ +P2,m-g, m^^, m^J cos' ^5 
+ — T^QhuQhAA-, ^T-^ Co(— Pi, Pi + P2, "^0, 1^1 .rnr) sin' 6*^ 

StT^^ '^^^(s-M^+iM^r^ i-l,i-l FI, b 



^9h22ghAAj^^^j^^^^—^^^Co ( -Pi , Pi + P2 , , mg^ , mg^ ) sin 9^ cos ^5 
^94 = -^^ftiiffftAA ^^ _ ^ iMhVh ^^^~^^'^^ +P2,mg, fflfc,, mgj sin cos 
+ ^gfei2gfe^A _ ^2^^^ iM^r^ ^"^"^''^' + P2, mg, m^^, m^J sin' 
- ^9hi2ghAA iMhVh ^^^'^^'^^ +P2,mg, mi,^,m-J cos' 0^ 

1 TTl ~ 

+ ^gh229hAA _ ^ iM;,rfe ^° ^ , Pi + P2 , , m^^ , mg^ ) sin 9^ cos ^5 
Xioi = XQi{ghjk ^ ffhAA ^ « = 3, 4; j, /c = 1, 2 

(B2) 
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APPENDIX C: HIGGS COUPLINGS IN MSSM 



9hbb 


giibb 


9 Abb 


gmi, sin a 
2Mw cosf) 


gmb cos a 
2Mw cosf) 


2Mw ^ 



TABLE I: The (j)bh vertex couplings. For (f> = h,H the Feynman rule is igifibb and for (p = A the 
Feynman rule is 'j^gAbb- 



9 


h 


H 


hh 


-32SF^°«2asin(/3 + a) 


-2^^[2sin2asin(/3 + Q) - cos(/3 + a) cos 2a] 


HH 


2c^l^ [2 sin 2a cos (/3 + a) + sin(/3 + a) cos 2a] 


-32^cos2acos(/3 + a) 


AA 


-2&^os2Psm{P + a) 


2^cos2/3cos(/3 + a) 



TABLE IL The 4'i(/)j<pj triple Higgs vertices. For (f) = h, H, A the Feynman rule is ^g^^(j,j(f,y 
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ghii 


2Mw cos p [sin i^i COS a + Absma)] + ^ sm a 
sin(/3 + a) [(1 ^^^3 ) cos^ 0^ + 2 tan^ 6*^/] 


gh22 


2Mw cos /3 [ sm 2^6 (/X COS « + ^fe sm «)] + ^^^^ ^ sm a 
2^ sin(/? + a) [(1 ^^^^3 ) sin2 6*^ + ^ tan^ Ow] 


9hl2 


2mTcos 13 COS 26lb(^ COS a + Ab sm a)] 
+ ^'^^ sin20fcSin(/? + a)(l ^^'^g^™) 


9HU 


2a4 cos /3 (sm 206) (// sm a cos a)] ^ cos a 
+^ cos(/? + a) [(1 '^'^l ) cos2 eb + ^ tan2 ^v^] 


9H22 


2Ml,cos/3Sin206(//sma A cos a) cos /3 cos a 
1 3*2'^cos(/3 1 a)[(l *'''^3^»)sin2 6'fe I ^ tan^ %] 


9H12 


ntf"^'' — 3 COS 20;, (/i sin a — Ah cos a) 

2Mvi/cos/3 "o\t^^ " 

^^^^ sin 20b cos(^ + q) (1 *^°3 ) 


9A12 


2Mrcos/3(Mcos/3 + Asm/3) 


gA21 


-5^12 



TABLE III: Higgs-sbottom-sbottom couplings. For (p = h,H the Feynman rule is ig^i i. and for 
(j) = A the Feynman rule is gj^i i - 
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